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Abstract 

A hierarchical, reversible mapping between levels of tree structured compu- 
tation, applicable for structuring the Quantum Computation algorithm for 
AAP-complete problem is presented. It is proven that confining the state of a 
quantum computer to a subspace of the available Hilbert space, where states 
are consistent with the problem constraints, can be done in polynomial time. 
The proposed mapping, together with the method of the state reduction can 
be potentially used for solving AAP-complete problems in polynomial time. 
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Motivation Expected running time scaling polynomially with input size, displayed by 
Shor's factoring algorithm [Q and, to even greater extent, the more recent results [0,^ give 
evidence that Quantum Mechanics offers a computation model more powerful than a classical 
one, yet still physically feasible. 

An application this power for simulating nondeterminism associated with ATP-complete 
problems had been discouraged by Bennett et al. in . However, Theorem 2 of 0] , which can 
be interpreted as a "no-go" proof for unstructured attempts to solve A/'P-complete problems 
in polynomial time, assumes particular quantum computation consisting of unitary evolution 
of simple input state followed by a single measurement of a state of a system. 

We will follow different approach, patterned after error correction, with repeated mea- 
surements forcing projections of a system state onto a subspace spanned on eigenvectors of 
particular eigenvalue of the Hermitian operator representing measurement. 

In remaining part of this paper, after choosing a known A/'P-complete problem as our 
object of study, we will present the deterministic, data independent mapping between the 
levels of a solution construction tree, thus providing a method for structuring quantum 
algorithms for ATP-complete problems. 

Sample Problem Let us consider a A/'P-complete combinatorial optimization problem 
known as a symmetric Traveling Salesman Problem (TSP) defined as follows: 
Given a weighted complete graph G = (V, E , w) , where \^ is a set of vertices Vk, k = l..n, 
E is a set of edges, Vfj,ffe G V {vi,Vk} G E and a function w : E ^ Z is symmetric: 
Wi^k = w{ei^k) = w{ek,i), 

find a Hamiltonian cycle t (a closed tour through the graph, visiting each vertex exactly 
once) with the minimal total weight: 

n-l 

W{t) = w{{Vt(^n),Vt(l)}) + w{{Vt(^k), Vtik+1)}) (1) 

k=l 

Instead of coding each cycle as a list of vertices, we will code it as a set of its edges, 
since in this case ordering information is coded locally, reducing number of bits modified at 
a cycle modification. We can represent each edge present in a cycle by a single set bit in 
a list of all n(n — l)/2 edges. The ordering of the list is such that the edges are grouped 
by the higher vertex number present in the edge and, within the group, by the number of 
the other vertex of the edge: (62,1,63^1,63^2,64,1, . . .) . In such a way, incomplete potential 
solutions i.e. cycles over first m vertices of a n- vertex problem have all bits above m(m+ 1)/2 
reset to 0. We need O^n"^) bits to code a n-vertex cycle, as compared with C(?T,log(n)) of 
more common coding as a list of vertices. Such a redundancy can be used for rudimentary 
error control, but the main motivation is robustness, universality and particular simplicity 
of unitary operator used for mapping between levels of a construction tree. 

From each cycle of m vertices we can construct m cycles connecting m + 1 vertices by 
breaking it at one of its m edges, and inserting the vertex m + 1 at the breaking point. We 
reset the relevant bit to zero (since the edge 6j,j we break at no longer exists in a resulting 
cycle ) and replace it by two bits representing edges which connect broken ends of the cycle 
of a previous level to the newly inserted vertex: Cm+ij, 6^+1,1- 

Reversible Mapping for TSP A quantum computation method of producing all the 
permutations of a given set of elements (all vertices Vi G V in a. case of TSP) has been 
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already provided in however, it requires resources scaling exponentially with number of 
elements in a set. 

It is interesting to note that latest results in Quantum Computation are not presented 
in Quantum Turing Machine formalism. This, in author opinion, reflects the growing belief 
that if we are ever able to build a quantum computer it will rather be in a form of uniform 
array of identical, locally interacting elements with quantum degrees of freedom than 
the counterpart of a classical computer. We will follow this approach. 

We start with the only tour connecting three vertices represented by a state |l2,il3,il3,20), 
where is used to denote that all the remaining qbits in all remaining groups are 0. In 
subsequent equations we will drop the indices for shorter notation, while preserving the 
same convention for the order of qbits. 

Let us denote /— th bit of a pure state \s > by s{l}. Clearly, the mapping of the previous 
paragraph is not unitary, but it can be made so by supplying an ancilla prepared in an equal 
superposition of all natural unit states u^il'} = with /, = l..m{m — l)/2. [^] 

Let us assume for a moment that our system is prepared in an equal superposition of all 
Hamiltonian cycles over first m vertices: 

Qbit set to 1 in ancilla {u]^{l} = 1) points to the edge in a tour of the input level which 
should be broken (if present i.e. if p^{l} = 1) by the insertion of a vertex of the next level. 
On output both these qbits are reset to 0, while setting two bits representing newly created 
edges p^,~^^{i} = P^i~^^{i'} = 1- All other bits of the input path remain unchanged, since the 
edges they represent are common for and p^i~^^- Partially expanded construction tree is 
depicted in Fig. |l], with levels labeled by the number of vertices. 

Thus defined Um is an identity operator except for the matrix elements: 

{uT,pT\0,Pr') = {0,p^^'\uT,pT) = 1, (3a) 
(O,p^+i|0,p™+0 = {u^,pTK,P^) = 0, (3b) 

for such k.k' that: 

where, if / = l..m(m — l)/2 is a position of a qbit representing the edge ei'^i", are the 
positions of qbits representing the edges em,i", em,i' and all other bits of p"^ are equal to 
respective bits of p^,~^^ and denotes that all the bits of ancilla are 0. 

Note that within an adopted convention there is one-to-one correspondence between 
indices of a given edge and a position of the relevant bit in a list, which can be determined 
in polynomial time: 

In order to find indices i, k of an edge, coded as an /-th bit, add the integers 1 + 2 + . . . 
until the sum is greater than or equal to I. Number of terms in this sum plus one equals i. 
Subtracting the last term from the sum and subtracting the result from / yields k. 
To complete the construction, we will use the following example: 
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Example 1. We want to construct the unitary transformation, mapping the only cycle 
over three vertices: |1110) to the superposition of three possible tours connecting four 
vertices using ancilla: 

|7i3) = ;^(|100) + |010) + |001)) (4) 

and acting with unitary transformation Um in a product of Hilbert spaces of our primary 
set of qbits Hp and ancilla Hu- 

i=(|100) + 1010) + |001))|1 11 0) ^ i=|000) X 

x(|0 11 110 0) + |1 01 101 0) + |1 10 001 0)) (5) 

We can construct this mapping in polynomial time for arbitrary m, factoring Um into 
m(m — 1)/2 operations [/^ by iterating over bits of ancilla and, for each 1 in turn, performing 
the four bit operation on affected qbits of p and a bit of ancilla as defined by relations 
of Eqs. ^. This construction for the above example is depicted in Fig. ^. 

Thus, in order to recursively implement this construction to produce the superposition 
of all cycles over n vertices we need the time C(n^). If we do not reuse the ancillae the space 
requirements will also be 0{n^). 

Measurement Note that the example describes the only case when the valid tour has 
all bits set, and so the ancilla has just right number of states to pass all the amplitude to 
outputs under operation of Us. 

In all subsequent levels some portion of the amplitude will remain in the input level. 
We can correct this by the measurement with projector = \0u){0u\ onto state |0„) of 
ancilla after using Um- Since all cycles of the output level are entangled with |0„) the perfect 
measurement does not destroy the superposition of the output tours, while reducing to 
the amplitude remaining in the input level tours. 

To establish the scaling behavior of the measurement part of the procedure let us consider 
in detail the apphcation of Um to the state of our system 

^ ^ m(m-l)/2 (m-l)!/2 

1=1 k=l 



m{m—l) {m—l)\ 
2 V 2 



■\J{m — l)m! 



Eio-0 + EK-,pr)| (6) 

c l,k 



where k is indexing all the states fulfilling Eqs. 0. They are m qbits set to 1 in each 

so they are m(m — l)!/2 = m!/2 such states — which can be expected, since it is the 
number of distinct cycles at the m + 1 level. The indices I run through the states where 
the set bit of ancilla corresponds to unset bit of . 

The measurement with Mm reduces the state of our primary register to: 
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compatible with the ancilla in a one dimensional subspace of Tiu defined by the eigenvector 
belonging to the eigenvalue 1 of the projector. The probability of finding the system in the 
desired state is: 



m{m — l)/2 m — 1 

If a measurement of any state of ancilla requires time T, the measurement with Mm will 
yield the result in average time t = T/Pm = (m — l)T/2. 

Thus, total measurement time for the construction of n- vertex cycles is (9(n^). To 
simplify the measurement we can use an auxiliary qbit la^) holding the result of a function 
AND{pm{l},u]^{l}) (represented by a square in Fig. |^) calculated for each As 
far as the state reduction is concerned, the measurement with is equivalent to the 

measurement with Mm- 

The reversibility of our construction has special meaning, due to the irreversible mea- 
surements Mm- application of to the final state \0,p"^~^^) results in each entangled 
with all uf^ representing set qbits of the cycle, but nevertheless can be used to push the am- 
plitude back to a single state in \p), while leaving the information of the cycles in (entangled) 
ancillae. 

The present construction is not time nor space optimal, but it is very simple — it is not 
using all the resources postulated for quantum computation so far and confines much of a 
calculation to isolated space, thus reducing errors. We need only real amplitudes, and our 
dynamics is truly digital: matrix elements of our unitary operators are only or 1. What is 
more, the polynomial amplitude errors will lead only to the polynomial increase of expected 
running time, and changing of phases between the states of a given level will not affect 
the computation at all. 

In principle we could have calculated the total weight for each path, but the probability of 
finding a minimum in single measurement will be vanishingly small for large n. We could have 
used ancillae better matching the requirements of unitarity and, instead of the postulated 
measurement, the measurement with a sum of elementary projectors onto all eigenvectors of 
associated operator (repeating Um with new ancilla n case that the resulting state of ancilla 
is not |0„), much like in quantum error correction). 

However, the application of the described mapping and postulated method of state re- 
duction will be very useful for polynomial time algorithms for ATP-complete problems, to 
be presented in a follow-up paper. 

The author wishes to thank N. Margolus for guidance into a field of reversible computa- 
tion and his hospitality at MIT, where this work originated and to gratefully acknowledge 
a valuable discussion with W. Zurek. 
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FIGURES 
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FIG. 1. Partially expanded solution construction tree for TSP, with Hamiltonian cycles coded 
by set qbits representing their edges, where levels of a tree are labeled by the number of vertices 
m. Um is a sought mapping between levels. is used to denote that all other qbits of a state 
are 0. Ancillae, necessary at all levels for unitarity of Um are not shown. 
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FIG. 2. Factoring of Um into the bitwise operations in the case of m = 3. Fuh circles and 
crosses denote respectively change from to 1 and from 1 to according to Eqn.^. For each of 
thus defined 4-bit operations affected bits of (positions of circles) can be determined from the 
set bit of the ancilla ui (a position of a cross) in polynomial time. Dashed line and a squares denote 
a possible extention of Um and an associated auxiliary qbit used to faciliate the measurements, 
necessary at all other levels of the construction. 
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